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Abstract 

We review the integrable structure of the Dirichlet boundary problem in two dimensions. 
The solution to the Dirichlet boundary problem for simply-connected case is given through 
a quasiclassical tau-function, which satisfies the Hirota equations of the dispersionless Toda 
hierarchy, following from properties of the Dirichlet Green function. We also outline a 
possible generalization to the case of multiply-connected domains related to the multi- 
support solutions of matrix models. 
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1 Introduction: the Green function and the Hadamard formula 



Solving the Dirichlet boundary problem PP, one reconstructs a harmonic function in a bounded 
domain from its values on the boundary. In two dimensions, this is one of standard problems 
of complex analysis having close relations to string theory and matrix models. Remarkably, 
it possesses a hidden integrable structure 2 . It turns out that variation of a solution to the 
Dirichlet problem under variation of the domain is described by an infinite hierarchy of non- 
linear partial differential equations known (in the simply-connected case) as dispersionless Toda 
hierarchy. It is a particular example of the universal hierarchy of quasiclassical or Whitham 
equations introduced in E] . 

The quasiclassical tau-function (or its logarithm F) is the main new object associated with 
a family of domains in the plane. Any domain in the complex plane with sufficiently smooth 
boundary can be parametrized by its harmonic moments and the F-function is a function of 
the full infinite set of the moments. The first order derivatives of F are then moments of 
the complementary domain. This gives a formal solution to the inverse potential problem, 
considered for a simply-connected case in [SJ El • The second order derivatives are coefficients 
of the Taylor expansion of the Dirichlet Green function and therefore they solve the Dirich- 
let boundary problem. These coefficients are constrained by infinite number of universal (i.e. 
domain-independent) relations which, unified in a generating form, just constitute the disper- 
sionless Hirota equations. For the third order derivatives there is a nice "residue formula" 
which allows one to prove [Jj that F obeys the Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) 
equaitons [H]. 

Let us remind the formulation of the Dirichlet problem in planar domains. Let D c be a 
domain in the complex plane bounded by one or several non-intersecting curves. It will be 
convenient for us to realize the D c as a complement of another domain, D (which in general 
may have more than one connected components), and consider the Dirichlet problem in D c . 
The problem is to find a harmonic function u(z) in D c , such that it is continuous up to the 
boundary <9D C and equals a given function uo(£) on the boundary, and it can be uniquely solved 
in terms of the Dirichlet Green function G(z,£): 



where d n is the normal derivative on the boundary with respect to the second variable, and the 
normal vector n is directed inward D c , \d£\ := dl(£) is an infinitesimal element of the length of 
the boundary <9D C . 

The Dirichlet Green function is uniquely determined by the following properties pQ: 

(Gl) The function G(z, z') is symmetric and harmonic everywhere in D c (including oo if D c 3 
oo) in both arguments except z = z' where G(z, z') = log \z — z'\ + . . . as z — > z'; 

(G2) G(z, z') = if any of the variables z, z' belongs to the boundary. 

Note that the definition implies that G(z,z') < inside D c . In particular, d n G(z,^) is strictly 
negative for all £ € (9D C . 

If D c is simply-connected (the boundary has only one component), the Dirichlet problem is 
equivalent to finding a bijective conformal map from D c onto the unit disk or any other reference 
domain (where the Green function is known explicitly) which exists by virtue of the Riemann 
mapping theorem. Let w(z) be such a bijective conformal map of D c onto the complement to 




(1.1) 
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Figure 1: A "pictorial" derivation of the Hadamard formula. We consider a small deforma- 
tion of the domain, with the new boundary being depicted by the dashed line. According 
to (G2) the Dirichlet Green function vanishes G(z,£) = if £ belongs to the old bound- 
ary. Then the variation 5G(z, £) simply equals to the new value, i.e. in the leading order 
8G(z,£) = —8n(£)d n G(z, . Now notice that 6G(z,£) is a harmonic function (the logarithmic 
singularity cancels since it is the same for both old and new functions) with the boundary value 
—5n(£)d n G(z,£). Applying ifl.l)) one obtains i)1.4|) . 



the unit disk, then 



G(z,z') = log 



w(z) — w(z') 



(1.2) 



w(z)w(z') — 1 

where bar means complex conjugation. It is this formula which allows one to derive the Hirota 
equations for the tau-function of the Dirichlet problem in the most economic and transparent 
way 0. Indeed, the Green function is shown to admit a representation through the logarithm 
of the tau-function of the form 



G(z,z') = log 



+ ±V(z)V(z')F 



;i.3) 



where V(z) (see (|2.5j) below) is certain differential operator with constant coefficients (depending 
only on the point z as a parameter) in the space of harmonic moments. Taking into account that 
G(z, oo) = — log |u;(,z)|, one excludes the Green function from these relations thus obtaining a 
closed system of equations for F only. 

Our main tool to derive (jl.3j) is the Hadamard variational formula 9 which gives variation 
of the Dirichlet Green function under small deformations of the domain in terms of the Green 
function itself: 

5G(z,z') = ^<f d n G(z,0d n G(z',05n(0\dH\. (1.4) 

ZTT JdD c 

Here 6n(£) is the normal displaycement (with sign) of the boundary under the deformation, 
counted along the normal vector at the boundary point £. It was shown in [2] that this remark- 
able formula reflects all integrable properties of the Dirichlet problem. An extremely simple 
"pictorial" derivation of the Hadamard formula is presented in fig. ^ Looking at the figure and 
applying (jl.ljl . one immediately gets (|1.4|) . 
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2 The Dirichlet problem for simply-connected domains and dis- 
persionless Hirota equations 



Let D be a connected domain in the complex plane bounded by a simple analytic curve. We 
consider the exterior Dirichlet problem in D c = C \ D which is the complement of D in the 
whole (extended) complex plane. Without loss of generality, we assume that D is compact and 
contains the point z = 0. Then D c is an unbounded simply-connected domain containing oo. 



2.1 Harmonic moments and elementary deformations 

To characterize the shape of the domain D c we consider its moments with respect to a complete 
basis of harmonic functions. The simplest basis is {z~ k }, {z~ k } (k > 1) and the constant 
function. Let t}~ be the harmonic moments 

t k =-—[ z- k d 2 z, fc = l,2,... (2.1) 



nk Jd 

and ik be the complex conjugated moments. The Stokes formula represents them as contour 
integrals t k = § dD z~ k zdz, providing, in particular, a regularization of possibly divergent 
integrals (|2.1|) . The moment of constant function is infinite but its variation is always finite 
and opposite to the variation of the complimentary domain D. Let to be the area (divided by 
vr) of D: 

to = - [ d 2 z (2.2) 



vr jd 

The harmonic moments of D c are coefficients of the Taylor expansion of the potential 

$( z ) = — / log \z - z'\d 2 z' (2.3) 
vr Jo 

induced by the domain D filled by two-dimensional Coulomb charges with the uniform density 
— 1. Clearly, d z dg$(z) = — 1 if z G D and vanishes otherwise, so around the origin (recall that 
D 3 0) the potential is — \z\ 2 plus a harmonic function: 

<&(*) - HO) = -\z\ 2 + ]T (t k z k + i k z k ) (2.4) 

k>\ 

A simple calculation shows that t k are just given by (|2.1j) . 

The basic fact of the theory of deformations of closed analytic curves is that the (in general 
complex) moments {th,i k } = {t±t} supplemented by the real variable to form a set of local 
coordinates in the "moduli space" of smooth closed curves. This means that under any small 
deformation of the domain the set t = {to, t± k } is subject to a small change and vice versa. For 
more details, see The differential operators 

V(z) = d t0 + £ ( Z -^d tk + ^dz k ) (2.5) 



k>l 



span the complexified tangent space to the space of curves. The operator V(z) has a clear 
geometrical meaning. To clarify it, we introduce the notion of elementary deformation. 

Fix a point z G D c and consider a special infinitesimal deformation of the domain such that 
the normal displaycement of the boundary is proportional to the gradient of the Green function 
at the boundary point (fig. EJ: 

Sn(0 = ~d n G(z,0 (2.6) 
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Figure 2: The elementary deformation with the base point z. 



For any sufficiently smooth initial boundary this deformation is well-defined as e —* 0. We call 
infinitesimal deformations from this family, parametrized by z £ D c , elementary deformations. 
The point z is refered to as the base point of the deformation. Note that since d n G < (see the 
remark after the definition of the Green function in the Introduction), Sn for the elementary 
deformations is either strictly positive or strictly negative depending of the sign of the e. 

Let 5 Z be variation of any quantity under the elementary deformation with the base point 
z. It is easy to see that S z to = e, 8 z t^ = ez~ k /k. Indeed, 



by virtue of the Dirichlet formula 

Let X = X(t) be any functional of our domain that depends on the harmonic moments 
only (in what follows we are going to consider only such functionals). The variation 5 Z X in the 
leading order in e is then given by 



The right hand side suggests that for functionals X such that the series V(z)X converges 
everywhere in D c up to the boundary, 5 Z X is a harmonic function of the base point z. 

Note that in [2] we used the "bump" deformation and continued it harmonically to D c . So 
it was the elementary deformation (J2Z3) S z oc § \d£\d n G(z,£)5 bump (£,) that was really used. The 
"bump" deformation should be understood as a (carefully taken) limit of S z when the point z 
tends to the boundary. 

2.2 The Hadamard formula as integrability condition 

Variation of the Green function under small deformations of the domain is known due to 
Hadamard, see eq. (|1.4j) . To find how the Green function changes under small variations of the 
harmonic moments, we fix three points a, b, c E C \ D and compute 5 c G(a, b) by means of the 
Hadamard formula (|1.4j) . Using 1)2.8 J) , one can identify the result with the action of the vector 
field V(c) on the Green function: 




<fc k Sn(0\d£\ = -^jc k d n G{z^)\di\ = - k z- k 



(2.7) 



5 z X = Y J ^ r 5 z t k = eV{z)X 



(2.8) 




(2.9) 
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Remarkably, the r.h.s. of ()2.9|) is symmetric in all three arguments: 

V(a)G{b, c) = V(6)G(c, a) = V(c)G(a, 6) (2.10) 

This is the key relation, which allows one to represent the Dirichlet problem as an integrable 
hierarchy of non- linear differential equations [2] ■ This relation is the integrability condition of 
the hierarchy. 

It follows from (|2.10|) (see [2] for details) that there exists a function F = F(t) such that 

+ l -V(z)V{z')F (2.11) 

The function F is (logarithm of) the tau- function of the integrable hierarchy. In ^3] it was 
called the tau-function of the (real analytic) curves. Existence of such a representation of the 
Green function was first conjectured by Takhtajan. This formula was first obtained in ^3] (see 
also 12 for a detailed proof and discussion). 



G(z,z') = \og 



2.3 Dispersionless Hirota equations for F 



Combining 1)2.11(1 and (jl.2j) . we obtain the relation 



log 



w(z) — w(z') 



w(z)w(z') — 1 



2 


1 1 


= log 






z z' 



+ V(z)V(z r )F 



(2.12) 



which implies an infinite hierarchy of differential equations on the function F. It is convenient 
to normalize the conformal map w(z) by the conditions that w(oo) = oo and d z w(oo) is real, 
so that 

w(z) = - + 0(l) as^oo (2.13) 
r 

where the real number r = lim^oo dz/dw{z) is called the (external) conformal radius of the do- 
main D (equivalently, it can be defined through the Green function as logr = \\m z ^ 00 {G{z^ oo) + 
log |z|), see ^U). Then, tending z' — > oo in (|2.12|) . one gets 



log |^(-z)| 2 = log |z| 2 — dt V(z)F 
The limit z — > oo of this equality yields a simple formula for the conformal radius: 

logr 2 = a 2 Q F 



(2.14) 



(2.15) 



Let us now separate holomorphic and antiholomorphic parts of these equations. To do that 
it is convenient to introduce holomorphic and antiholomorphic parts of the operator V(z) (|2.5j) : 



DM 



D(z) 



k>l 



Rewrite (|2,12jl in the form 

w{z) — w(z') 



log 



w(z)w(z') — 1 / 



z z' 



(2.16) 



fe>i 



-d t0 +D(z))V(z')F 



log 



w(z) — w(z') 
w(z')w(z) — 1 



+ log 



1 1 



z z 



+ [^d t0 +D(z))V(z')F 
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The l.h.s. is a holomorphic function of z while the r.h.s. is antiholomorphic. Therefore, both 
are equal to a z-independent term which can be found from the limit z — > oo. As a result, we 
obtain the equation 

^;^w)-*K) +DWVM * (2 ' 17> 

which, as z' — > oo, turns into the formula for the conformal map 



log«;(2) = log* - -5 t 2 F - a t()J D(z)F (2.18) 



(here we used (|2.15[> ). Proceeding in a similar way, one can rearrange (|2.17|) in order to write 
it separately for holomorphic and antiholomorphic parts in z': 



*>(*)-">(*') = _i_ <F + D(z)D(z , )F (2 . 19) 



- log (1 - ' ) = D(z)D(z>)F (2.20) 

\ W[Z)W{Z ) ) 

Writing down eas. (|2.19|) for the pairs of points (a, b), (b,c) and (c, a) and summing up the 
exponentials of the both sides of each equation one arrives at the relation 

(O _ b)e D(a)D(b)F + {h _ c)e D(b)D(c)F + ((J _ a)e ^( C )D(a)F = q (2 21) 

which is the dispersionless Hirota equation (for the KP part of the two-dimensional Toda lattice 
hierarchy) written in the symmetric form. This equation can be regarded as a very degenerate 
case of the trisecant Fay identity. It encodes the algebraic relations between the second order 
derivatives of the function F. As c — > oo, we get these relations in a more explicit but less 
symmetric form: 

1 _ e D(a)D(b)F = D(a) ~ DiP) F 

a — b 

which makes it clear that the totality of second derivatives i 7 ^- := d^dt-F are expressed through 
the derivatives with one of the indices equal to unity. 

More general equations of the dispersionless Toda hierarchy obtained in a similar way by 
combining eqs. Q2.18JI . ()2.19|) and (|2.2U|) include derivatives w.r.t. to and 

( a _ b)e D(a)D(b)F = ae -d t0 D(a)F _ be -d to D(b)F {2 2?>) 

1 _ e -D(z)D(z)F = J_ e d tQ V(z)F ( 2 _ 24 ) 

zz 

These equations allow one to express the second derivatives dt m dt n F, dt m di n F with m,n > 1 
through the derivatives dt dt k F, dt df k F. In particular, the dispersionless Toda equation, 

dtidhF = ( ? 2 t F (2.25) 

which follows from Q2.24j) as z — > oo, expresses d^d^F through df Q F. 

For a comprehensive exposition of Hirota equations for dispersionless KP and Toda hierar- 
chies we refer the reader to [151 116| . 
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2.4 Integral representation of the tau-function 



Eq. (j2.ll.jl allows one to obtain a representation of the tau-function as a double integral over the 
domain D. Set &(z) := V(z)F. One is able to determine this function via its variation under 
the elementary deformation: 

6 a $(z) = -2elog a' 1 - z' 1 + 2eG(a,z) (2.26) 

which is read from eq. (j2.11j) by virtue of (j2,8j) . This allows one to identify $ with the "modified 
potential" <&(z) = &(z) — <I>(0) + ^ log \z\ 2 , where $ is given by (j2.3j) . Thus we can write 



V(z)F = $(z) =-\\ lo § \ z ~ l ~ C VC = v + 2ReY, T (2.27 ) 



k>0 



V ±„-k 
k 



The last equality is to be understood as the Taylor expansion around infinity. The coefficients 
Vk are moments of the interior domain (the "dual" harmonic moments) defined as 

v k = ~ [ z k d 2 z (k>0), v = -$(0) = - / log \z\d 2 z (2.28) 

From (j2.27j) it is clear that 

v k = d tk F, k>0 (2.29) 
In a similar manner, one obtains the integral representation of the tau-function 

F = ~ f [ \og\z- 1 -C X \d 2 zd 2 C (2.30) 

or 

F = — f $(z)d 2 z (2.31) 



2vr 

These formulas remain intact in the multiply-connected case (see below). 



3 Towards the multiply-connected case and generalized Hirota 
equations 

Now we are going to explain how the above picture can be generalized to the multiply-connected 
case. The details can be found in j!7j . 

Let D a , a = 0, 1, . . . , g, be a collection of g + 1 non-intersecting bounded connected domains 
in the complex plane with smooth boundaries dD a . Set D = U^, = qDq,, so that the complement 
D c = C \ D becomes a multiply-connected unbounded domain in the complex plane (see fig. EI), 
b a being the boundary curves. 

It is customary to associate with a planar multiply-connected domain its Schottky double, a 
compact Riemann surface without boundary endowed with an antiholomorpic involution, the 
boundary of the initial domain being the set of fixed points of the involution. The Schottky 
double of the multiply-connected domain D c can be thought of as two copies of D c ("upper" 
and "lower" sheets of the double) glued along the boundaries U 9 a=0 b a = dD c , with points at 
infinity added (oo and 6b). In this set-up the holomorphic coordinate on the upper sheet is 
z inherited from D c , while the holomorphic coordinate on the other sheet is z. The Schottky 
double of D c with two infinities added is a compact Riemann surface of genus g = #{D Q } — 1. 
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Figure 3: A multiply-connected domain D c = C \ D for g = 3. The domain D = U«=o 
consists of g + 1 = 4 disconnected parts D a with the boundaries b a . To define the complete 
set of harmonic moments, we also need the auxiliary points z a € D a which should be always 
located inside the corresponding domains. 

On the double, one may choose a canonical basis of cycles. The b-cycles are just boundaries 
of the holes b a for a = 1, ... ,g. Note that regarded as the oriented boundaries of D c (not D) 
they have the clockwise orientation. The a a -cycle connects the a-th hole with the 0-th one. To 
be more precise, fix points £ a on the boundaries, then the a a cycle starts from £o> goes to £ Q 
on the "upper" (holomorphic) sheet of the double and goes back the same way on the "lower" 
sheet, where the holomorphic coordinate is z. 

3.1 Tau-f unction for algebraic domains 

Comparing to the simply-connected case, nothing is changed in posing the standard Dirichlet 
problem. The definition of the Green function and the formula (jl.lj) for the solution of the 
Dirichlet problem through the Green function are the same too. A difference is in the nature of 
harmonic functions. Any harmonic function is the real part of an analytic function but in the 
multiply-connected case these analytic funstions are not necessarily single-valued (only their 
real parts have to be single- valued) . 

One may still characterize the shape of a multiply-connected domain by harmonic moments. 
However, the set of linearly independent harmonic functions should be extended. The complete 
basis of harmonic functions in the plane with holes is described in [T7J. 

Here we shall only say a few words about the case which requires the minimal number of 
additional parameters and minimal modifications of the theory. This is the case of algebraic 
domains (in the sense of 1 1 j ) . or quadrature domains |181 119j . where, roughly speaking, the 
space of independent harmonic moments is finite-dimensional. For example, one may keep in 
mind the class of domains with only finite number of non- vanishing moments. It is this class 
which is directly related to multi-support solutions of matrix models with polynomial poten- 
tials. Boundaries of such multiply-connected domains can be explicitly described by algebraic 
equations |2flj . 

In this case it is enough to incorporate moments with respect to g additional harmonic 
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functions of the form 



u a (z) = log 



Z a 
Z 



a = l,...,g 



where z a € D a are some marked points, one in each hole (see fig- Ej) - Without loss of generality, 
it is convenient to put zq = 0. The "periods" of these functions are: §^ d n up(z)\dz\ = ^Tr5 a p. 
The independent parameters for algebraic domains are: 

1 r , 9 Area(D) 
t = - d 2 z- 



vr Jo vr 



t k = ~r f z- k d 2 z, k > 1 (3.1) 



Z 2 

log 1 

(tfc are complex numbers while to an d 4> a are real). Instead of <f) a it is more convenient to use 



- / log 



d 2 z , a = 1, . . . , g 



n a = 4> a - 2 Re tkZa (3.2) 
fc>o 

which does not depend on the choice of z a 's. Note that in the case of a finite number of 
non-vanishing moments the sum is always well-defined. 

Using the Hadamard formula, one again derives the "exchahge relations" l)2.1fl|) which imply 
the existence of the tau-function and the fundamental relation (|2.11|) . They have the same form 
as in the simply-connected case. It can be shown that taking derivatives of the tau-function 
with respect to the additional variables Uq,, one obtains the harmonic measures of boundary 
components and the period matrix. 

The harmonic measure uo a (z) of the boundary component b a is the harmonic function in D c 
such that it is equal to 1 on b a and. vaiiish.es on the other boundary curves. From the general 
formula f|l.ljl we conclude that 

u a (z) = -—I d n G(z,()\d(\, a = l,...,g (3.3) 
2vr Jb a 

Being harmonic, io a can be represented as the real part of a holomorphic function: 



u) a (z)=W a (z) + W a (z) 

where W a (z) are holomorphic multivalued functions in D c . The differentials dW a are holo- 
morphic in D c and purely imaginary on all boundary contours. So they can be extended 
holomorphically to the lower sheet of the Schottky double as —dW a {z). In fact this is the 
canonically normalized basis of holomorphic differentials on the double. Indeed, according to 
the definitions, 

dW p = 2Re / " dWpiz) = up(S a )-up(to) = 5 a(3 

Then the matrix of b-periods of these differentials reads 

T Ql 3 = f dWp = - % - f d n ujpdl = iirVL af3 (3.4) 

The period matrix T a p is purely imaginary non-degenerate matrix with positively definite imag- 
inary part. In addition to (|2.11l) . the following relations hold: 

io a (z) = -d a V{z)F (3.5) 
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and 

T Q = 2md a dpF (3.6) 

where d a := d/dU a . 

In the multiply-connected case, the suitable analog of the conformal map w{z) (or rather of 
log w(z)) is the embedding of D c into the g-dimensional complex torus Jac, the Jacobi variety 
of the Schottky double. This embedding is given, up to an overall shift in Jac, by the Abel 
map z i — > W(z) := (W±(z), . . . , W g (z)) where 

W Q (z) = f Z dW a (3.7) 
J to 

is the holomorphic part of the harmonic measure oj a . By virtue of (|3.5|h the Abel map is 
represented through the second order derivatives of the function F: 

W a {z) - W a (oo) = [ Z dW a = -d a D(z)F (3.8) 

J oo 

2 Re W a (oo) = co a (oo) = -d to d a F (3.9) 
The last formula immediately follows from (|3.5|) . 



3.2 The Green function and the generalized Hirota equations. 

The Green function of the Dirichlet boundary problem in the multiply-connected case, can be 
written in terms of the prime form (see for the definition and properties) on the Schottky 
double (cf. Cd)): 

^ C )=log|^ C) 



,,. , (3-10) 

Here by £ we mean the (holomorphic) coordinate of the "mirror" point on the Schottky double, 
i.e. the "mirror" of £ under the antiholomorphic involution. Using (|3.1Uj) together with (|2.11|) . 
(|3.5|) and (|3.6j) one can obtain the following representations of the prime form in terms of the 
tau-function 

E(z,C) = (z- 1 - (-^e-hi^-DiOfF 

iE(z,z) = e-^ 2 ^ F 

generalizing ()2.19|) . 1)2.20(1 in the simply-connected case. 

This allows us to write the generalized Hirota equations for F in the multiply-connected 
case. They follow from the Fay identities and (|3.11|) . In analogy to the simply-conected 
case, any second order derivative of the function F w.r.t. tj. (and F^, is expressed through 
the derivatives {F a p} where a, (3 = 0, ... ,g together with {F ati \ and their complex conjugated. 
To be more precise, one can consider all second derivatives as functions of {F a p, F^} modulo 
certain relations on the latter discussed in \TJ\; sometimes on this "small phase space" more 
extra constraints arise, which can be written in the form similar to the Hirota or WDVV 
equations [22]. 

For the detailed discussion of the generalized Hirota relations the reader is addressed to 
|17j . Here we just give the simplest example of such relations, an analog of the dispesrionless 
Toda equation (|2.25j) for the tau-function. It reads 

W = ^(~)+gff(~)- g O e ^ + ± ( lo g e(Z)U,(d a O tl F)( dp d h F) (3.12) 

6 W a,/3=l 
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Here is the Riemann theta-function with the period matrix T a p and 

(log0(Z)), a/3 := d 2 log9(Z)/dZ a dZ l3 



The equation holds for any vector-valued parameter Z G Jac. It is important to note that the 
theta-functions are expressed through the second order derivatives of F, so (|3.12J) is indeed a 
partial differential equation for F. For example, 



4 Conclusion 

In these notes we have reviewed the integrable structure of the Dirichlet boundary problem. 
We have presented the simplest known to us proof of the Hadamard variational formula and 
derivation of the dispersionless Hirota equations for the simply-connected case. 

We have also demonstrated how this approach can be generalized to the case of multiply- 
connected domains. The main ingredients remain intact, but the conformal map to the reference 
domain should be substituted by the Abel map into Jacobian of the Schottky double of the 
multiply-connected domain. Then one can write the generalization of the Hirota equations 
using the Fay identities, a particular case of which leads to generalization of the dispersionless 
Toda equation. 

Here we have only briefly commented on the properties of the quasiclassical tau-function 
of the multiply-connected solution. A detailed discussion of this issue and many related prob- 
lems, including conformal maps in the multiply-connected case, duality transformations on the 
Schottky double, relation to the multi-support solutions of the matrix models etc, can be found 
inHZI. 
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